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It is argued that the presence of a nonanalytic term in the effective potential of the Ginzburg-Landau model is 
immaterial as far as the order of the superconductor-normal phase transition is concerned. To achieve agreement 
with the renormalization group, the effective potential has to be extended to include derivative terms, which 
provide the theory with a low momentum scale which can be varied to probe the (possible) fixed point. 



I. PRELUDE 



This paper is dedicated to Professor Kleinert on the occasion of his 60th birthday with special thanks for the exciting years I 
had the good fortune to be his Wissenschaftlicher Assistent at the Freie Universitat Berlin. It is an implementation of one of the 
many neat ideas he shared with me over the years. This particular one he told me about four years ago in the late summer of 
^ ! 1996. 
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It is a commonly held belief that the order of an equilibrium phase transition can be inferred from the form of the effective 
potential at criticality. The presence of, for example, a cubic term in the order parameter is taken as signalling a discontinuous 
transition, while the absence of such a term is taken as indicating-.a continuous transition. At the mean-field level, where 
fluctuations in the order parameter are ignored, this is certainly truetl. However, the question about the order of a transition is 
settled in the full theory, and fluctuations may well change the mean-field result. 
' In a renormalization-group (RG) approach, a continuous phase transition is associated with an infrared-stable fixed point in the 
, space of coupling parameters characterizing the theory. Sufficiently close to the transition, such a fixed point acts as an attractor 
to which the couplings flow when one passes to larger length scales by integrating out field fluctuations of smaller length scales. 
Tj" \ When no infrared-stable fixed point is detected in this process, the transition is discontinuous. 

It is important to note that a fixed point is probed by changing a scale. The effective potential evaluated at criticality in itself 
is inadequate for this as it lacks a scale. In this contribution we show that a scale can be introduced by extending the effective 
potential to include derivative terms. The resulting effective free energy contains the same information as that obtained in RG. 
Along the way we are able to implement Kleinert's idea of finding a fixed point without calculating flow functions first. 
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II. CUBIC TERM 



To be specific we consider the superconductor-normal phase transition described by the O(n) Ginzburg-Landau model which 
has been one of Kleinert's research topics for many years, and to which a large part of the first volume of his textbookcl Gauge 
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q Fields in Condensed Matter is devoted. The model is specified by the free energy density (in the notation of statistical physics) 

E = \{d» - zeA^ 2 + m 2 \cp\ 2 + A|0| 4 + \f% + ^{d.Arf, (1) 
^ ■ with a complex order parameter <fi having an even number n of real field components: 
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A conventional superconductor with Cooper pairing in the s-channel corresponds to n = 2. The parameters e and m are the 
electric charge and mass of the field, while A is the coupling constant characterizing the 4th-order interaction term. We included 
a gauge-fixing term with parameter a, and F^ v — d^Au — d^A^ is the (magnetic) field strength. For convenience we will work 
in the gauge d^A^ = 0, which is implemented by taking the limit a — > 0. The mass term depends on the temperature T, and 
changes sign at the critical temperature T c . In the Ginzburg-Landau model, m 2 = £,q 2 {T/T c — 1), where £o is the length scale 
of fluctuations in the amplitude of the order parameter. 

As a .first step to investigate the effect of fluctuations on the mean-field picture, we consider, following Halperin, Lubensky, 
and MaQ, those in the vector field A^. Since (at least in the normal phase) the corresponding mode is gapless, it can have an 
important impact on the infrared behavior of the theory. The functional integral over A^ is a simple Gaussian and gives as 
contribution to the effective potential energy density in d dimensions 

Vef f = ^/^ln(fc 2 + 2e 2 |$| 2 ), 0) 
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FIG. 1: Sketch of the 1-loop effective potential. 



assuming that the order parameter is a nonfiuctuating background field denoted by $. In deriving this, it is used that the 
combination 

iV(fc) = <W-^, (4) 

which appears in intermediate steps, is a projection operator satisfying P 2 = P, and that its trace gives the number of transverse 
components: tr P = d — 1. The momentum integral is easily carried out using dimensional regularization, with the result 

Vcff ^ (2^V r(1 " d/2)e ^ l<i ' (5) 

where T(a;)_is the Gamma function. For d = 3, this leads to the cubic term in the effective potential mentioned in the 
introductionE The 1-loop contribution is to be added to the mean-field potential Vo = m|$| 2 + A|$| 4 , and is often taken as 
indicating a discontinuous phase transition at a — what KleinerO likes to call — precocious temperature T\ above the mean-field 
critical temperature T c , 

where the mass term is still positive, and the order parameter jumps from zero to the finite value |$| 2 = (l/187r 2 )e 6 /A 2 (see 
FigJ). 

Since this result is obtained in perturbation theory, both e and A are assumed to be small. This still leaves the ratio of the 
two, or the so-called Ginzburg-Landau parameter Kq L = e 2 /A undetermined. This parameter separates type-II (kql > l/\/2) 
superconductors, which have a Meissner phase where an applied magnetic field can penetrate the sample in the form of quantized 
flux tubes, and type-I (kql < 1/V2) superconductors, for which these flux tubes become unstable. Neglecting fluctuations in 
the order parameter, as is done here, is valid in the type-I regime, where kql is small. (In the opposite limit, of deep type-II 
superconductors where kql is large, fluctuations in the vector field.can be neglected instead.) This then led to the conclusion 
that type-I superconductors undergo a discontinuous phase transitionH. It should, however, be noted that fluctuations in the order 
parameter produce at criticality also a cubic term in the effective potential. According to this argument, this would imply that 
also type-II superconductors should undergo such a transition. 

As the result discussed in this section is independent of the number of field components, it should be valid for any n, including 
large numbers. This opens the possibility to check it using al/n expansion. 



m. l/n EXPANSION 

The l/n expansion can be applied when the number n of field components is large, so that its inverse provides the theory with 
a small parameter. Contributions are then ordered not according to the number of loops, as in the loop expansion, but according 
to the number of factors l/n. The leading contribution in l/n due to fluctuations in the vector field is obtained by dressing its 
correlation function with arbitrary many bubble insertions, and summing the entire set of Feynman diagramsEl The resulting 
series is a simple geometrical one, which leads to the following change in the correlation function: 

Pfiu(k) _^ P[iv(k) 
k 2 k 2 — bane 2 k d ~ 2 ' 



3 



where bd — Cd/{d — 1), with Cd the 1-loop integral 

d d k 



(27r) d k 2 {k + q)< 



q 2 = l 



r(2-d/2)r 2 (rf/2-i) 

(47r) d / 2 r(d - 2) ' ' 



Below d = 4, the infrared behavior has become less singular as a result of including the bubble insertions. Instead of the effective 
potential (|J), we now obtain 

V cff = ^| ^ln(fc 2 + & d e 2 fc d - 2 + 2e 2 |<i>| 2 ). (9) 

For large fields, |$| > (b d e d ~ 2 ) 1/{ ' i - d \ the main contribution to the integral comes from large momenta, where k 2 +bde 2 k 4 « 
k 2 , so that in this limit we recover the result (Q) obtained in the loop expansion: V c ff oc In the limit of small fields, on the 
other hand, the main contribution comes from small momenta, where k 2 + bde 2 k d ~ A ss bde 2 k d ~ 4 . We then find V c ff oc |$| 2d 
instead, implying that in d = 3, the l/n expansion no longer gives a cubic term in the effective potential. [This argument 
captures only those terms in the Taylor expansion of the logarithm in the effective potential (^|) which diverge in the infrared. 
These are the important ones for our purposes as they can produce nonanalytic behavior. The first terms in the expansion may 
be infrared finite, depending on the dimensionality d, and have to be treated separately. But they always lead to analytic terms 
and are therefore of no concern to us here.] 

The l/n expansion calls into question the validity of the result obtained in the loop expansion as it corresponds to large field 
values and it is not clear whether this is still in the realm of perturbation theory. 

In the next section, we show that by extending the calculation of the effective potential to that of the effective free energy, 
which includes derivative terms, matters become consistent and in agreement with RG. 

IV. EFFECTIVE FREE ENERGY 

When computing the effective free energy, not only fluctuations in the vector field are to be considered, but also those in the 
order parameter. To this end, we set 

= $0) + $(x), (10) 

with $(x) a nonfluctuating background field, and integrate out the fluctuating field <fi. The main difference with the previous 
calculation of the effective potential is that the background field is no longer assumed to be constant, but can vary in space. We 
therefore also have to analyze the derivative terms in the effective theory. 

Because of gauge invariance, the results will depend solely on the absolute value |$|. Without loss of generality we can 
therefore assume that the background field $ has only one nonzero component, say the first one, $ = (v,0,--- ,0)/V2. For 
convenience we will work at criticality and set the mass parameter m to zero. 

A straightforward calculation yields in the gauge d^A^ = the 1-loop effective free energy 



F, n [r\ -- -Trln (l + 3A-^ 2 j + ^y^Trln (l + X^v 2 ) (II) 



•-Trln 
2 
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where we ignored an irrelevant constant, and the trace Tr denotes the trace tr over discrete indices as well as the integration over 
momentum and space. More precisely, 

Tr\n[l + K{x,p)]=tx J d d x J m [i + K ( X)P )] (12) 

with = —idfi the derivative operating on everything to its right. Since the background field is space-dependent, the integrals 
in Eq. ([[T]) cannot be evaluated in closed form, but only in a derivative expansioro 

The first step in this scheme is to expand the logarithm in a Taylor series. Each term of the series contains powers of the 
derivative which in the second step are shifted to the left using the identity 

f{x)Pix9{x) = - id ll )f(x)g(x), (13) 

where f(x) and g(x) are arbitrary functions. The symbol <9 M in the last term denotes the derivative which, in contrast to p^, 
operates only on the next object to its right. The next step is to repeatedly integrate by parts until all the p^ derivatives operate 
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to the left. They then simply produce factors of fc M as only the function exp( 
arbitrary function h(k): 



-ik 



appears at the left. In shorthand, for an 



d d k 



-ik-x 



h(p)f(x) = 



d d k 



-ik-x 



■ X h(-P)f(x) = 



d d k 



e- lk - x h(k)f(x), 



(14) 



ignoring total derivatives. In this way, all occurrences of the derivative operator are replaced with a mere integration variable 
The function exp(ik ■ x) at the right can now be moved to the left where it is annihilated by the function cxp(— ik ■ x). The 
momentum integration can then, in principle, be performed and the effective free energy takes the form of a space integral over 
a local density F e g — J d d x£ c g. 

Applied to the formal expression (|TT|), the derivative expansion yields for the quadratic terms in the expansion of the logarithms 
there, the following terms in the effective free energy density: 



feff(|*|) = 



-(n + 8)A 2 



d d k 1 



1 



) d k 2 (k - id) 2 
d d k i i 



(27r) d k 2 (k - id) 2 
d d k 1 1 



$|2|$|2 



P^{k)- 



(2n) d k 2 '{k-id) 



P vll {k - id)\Q\ 2 \®\ 



(15) 



Assuming that the order parameter carries a momentum k, we can replace all occurrences of the derivative — id^ with k^. 

It is important to note the absence of infrared divergences here since acts as an infrared cutoff. Without such a cutoff, as 
in the effective potential (^J), the individual terms in the Taylor expansion cannot be integrated because of infrared divergences, 
and the entire series has to be summed, which can lead to nonanalytic contributions. From the perspective of RG, the presence 
of the momentum scale is crucial as it allows the study of the fixed point by letting that scale approach zero. 

With the tree contribution £tree added, we obtain after carrying out the integrals over the loop momentum 



ftrce + £ c s = [1 - c (d - l)e 2 ] \d^\ 2 + {a - c [(n + 8)A 2 + d(d - l)e 4 ] } K 4 - d |$| 4 



(16) 



where A = \n d ~ A and e 2 = e 2 K d ~ 4 are the rescaled dimensionless coupling constants. As an aside, by evaluating the integrals 
in fixed dimension 2 < d < 4, and not-in an e expansion close to the upper critical dimension d = 4, we in effect implement 
Parisi's approach to critical phenomenalZl The factor in front of the kinetic term in Eq. (flq) amounts to a field renormalization. 

-1/2 

It can be absorbed by introducing the renormalized field $ r = $, with Z^, the field renormalization factor 

Z^ = l + c d {d-l)e 2 . 



(17) 



The effective free energy density then becomes 



4 

r| ) 



with A r the renormalized coupling 



1 _ 1 

A ~ A r ° d 



(„ + 8)-2(d-l)i + id(d-l)i 

A r 4 A r 



(18) 



(19) 



It was Kleinert's idea to put this equation in this form. The reason being that the critical point is approached by letting the 
momentum scale approach zero. Since the original coupling constant A is fixed, the left side tends to zero when n — > 0. The 
resulting quadratic equation then determines the value of A r at the critical point, provided we know the value of e 2 there. This 
procedure is equivalent to finding the root of the flow equation for A r in conventional RG0. 

To obtain the value of e 2 at criticality, we rescale the vector field — > A^/e and consider it instead of the order parameter 
to be the background field. The 1-loop effective free energy obtained after integrating out the scalar fields reads: 



-Tr In 



1 



1 

2 



4) 



(20) 



where we again ignored an irrelevant constant and used the gauge d^A^ = 0. The second term in the expansion of the logarithm 
yields the first nonzero contribution 



£ eS (A) = 



d d k kjjjsy - id v ) 
(2ir) d k 2 (k-id) 2 p 



(21) 
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If we assume that the vector field carries the same momentum as does the order parameter, we obtain after carrying out the 
momentum integral for the sum of the tree and the 1 -loop contribution 

£trcc + feff = -^d 2 A^, (22) 



with e 2 . the renormalized coupling constant, 

1 1 n Cd 



2d-l 



(23) 



Keeping e fixed and letting n — ► to approach the critical point, we see that the renormalized coupling tends to a constant value 

,9 2d-l 

e* 2 = . (24) 

n c d 

When this value is substituted in Eq. (|l9|) with the left side set to zero, the resulting quadratic equation in A r has two real 
solutions 

* = 2^)Z [n + A{d - 1?± ^ (25) 

provided that the determinant 



A d := v/n 2 - 4(d - 2)(d - l) 2 (d + ljn - 16(d - l) 3 (d + 1) (26) 

is real. For fixed dimensionality, this condition is satisfied only for a sufficient number of field components n. Specifically, the 
minimum number in d = 2, 3, 4 is n c (2) = 4^ « 6.9, n c (3) = 16(2 + VE) m 71.2, anJn c (4) = 12(15 + 4\/l5) « 365.9. 
The plus sign in Eq. ( |26"| ) corresponds to the infrared-stable fixed point, while the minus sign corresponds to the tricritical point, 
where the continuous phase transition changes to a discontinuous one. For a fixed number of field components n, this 1-loop 
result shows that continuous behavior is favored when d — > 2, while in the opposite limit, d — > 4, discontinuous behavior is 
favored. 



V. CONCLUSIONS 



For a conventional 3-dimensional superconductor corresponding to n = 2, the 1-loop result fails to give an infrared-stable 
fixed point. It should, however, be kept in mind that there are infinitely many loop diagrams, and the values of the coupling 
constants at criticality are not particular small to justify low-order perturbation theory. For example, e* 2 = 32/ro according to 
Eq. (|24j) with d = 3, which only for large n is small. Different approaches are therefore required to investigate the presence of a 
fixed point. |-_— .|— . 

One is the dual approacrfflEj, which is a formulation in terms of magnetic vortex loops. Monte Carlo simulations of the dual 
model led to the conclusion that 3-dimensional type-II superconductors undergo a continuous phase transition belonging to the 
XY universality class with an inverted temperature axisliil. Using this approach, KleinertO predicted the presence of a tricritical 
point at a value of the Ginzburg-Landau parameter K tri w 0.8/V2. ■— I 

Another approach within the framework of the Ginzburg-Landau model itself was put forward in Ref.U, where — in our 
language — the vector field was assumed to carry not the same momentum k m as the order parameter, but k^/x instead. The 
charged fixed points for n — 2 and d = 3 are then located at e* 2 = 16/a; which becomes small for large x, facilitating the 
existence of an infrared-stable fixed point. The free parameter x was determined by matching K% T i with Kleinert's estimate-pr 
Monte Carlo simulationsLj which gave K t ri ~ 0.42/V2. Form Eqs. (^4|) and < |23| > one obtains with the parameter x included^ 

K t 2 ri = (8 + x + Vx 2 + 16a; - 176)/40. (27) 

The resulting values for the critical exponents are consistent with the expected (inverted) XY universality class. 

Also resummation techniques applied to resultslij obtained to 2nd-order in the loop expansion predict the existence of an 
infrared-stable fixed point in the Ginzburg-Landau modeltZl. 

In conclusion, by extending the effective potential to include derivative terms, we achieved agreement with RG as they provide 
the theory with a low momentum scale which can be varied to probe the fixed point. The presence of a nonanalytic term in the 
effective potential at criticality was argued to be immaterial as far as the order of the phase transition is concerned. 
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